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1. Introduction. The Lebowitz inequality was first established for ferromagnetic Ising systems 
in [1]. It was later shown [2,3] that in these and related systems such an inequality implies a 
continuous transition. The inequality itself states that the (untruncated) four point correlation 
function is bounded above by the sum of the products of the two point functions paired in all 
possible ways. From this one can show that the derivative of the susceptibility 1 is bounded above 
by a constant times the square of the susceptibility. And from this one can show that if at some 
(3* the susceptibility, X, is infinite, then for any (3 < (3* we have X > [const.] [(3* — /3] _1 . This 
necessarily implies that the susceptibility has to diverge continuously at some point. 2 

In spin systems, the derivations of these sorts of inequalities usually follow the same course: 
Duplicate the system, rearrange variables, expand out all terms and hope everything comes out 
positive. (Although in [1], some FKG properties of the duplicated system were also exploited.) 
The model we consider, introduced in [6], is described by the Hamiltonian 



where a and r represent Ising variables and J and U are both positive. Although a proof 
along the above described lines may be possible, it would be extraordinarily complicated due 
to the presence of the negative couplings. (Not to mention the number of variables that would 
be required.) Our strategy - more akin to the derivation in [1] - involves two distinct steps. 
First: get a bound on the pure four-point function (involving only the a-variables) in terms 
of the various mixed four-point functions. Second: use the repulsive nature of the interaction 
between the two types of variables to bound the mixed four-point functions by products of the 
corresponding two-point functions. In short, no duplication. In point of fact, we will present 
these steps in the reverse order. In Section 2, we describe a random cluster representation for this 
system (introduced in [7,8]) that provides the bounds described in the second step. In Section 
3 we present a variant of the so-called loop expansion which gets us through the first step. We 
remark, without understanding, that both representations have their advantages and limitations. 
In particular, a unified proof using either representation alone seems to be completely hopeless. 
We remark without proof that this general strategy can be implemented in other systems as well, 
e.g. the usual 0(2) and 0(3) ferromagnets. However the system described by the Hamiltonian in 
Eq.(l.l) is the first new system for which a Lebowitz inequality can be proved by these methods. 

It is clear that the behavior of the system depends drastically on the magnitude of U. E.g. 
if U < J, the ground state is ferromagnetic and if U > J, it is not. We will work in the 
ferromagnetic domain which, it seems, is determined by the condition 



1 Throughout this work, we will use the word "susceptibility" to denote the sum of the untruncated 2-point 
function. For the systems under consideration, this coincides with the thermodynamic susceptibility in the single 
phase regime and it is the divergence of this object that we wish to study. 

2 This is easily seen on a formal level in infinite volume; a rigorous proof often involves finite-volume cutoffs. 
Such derivations are routine and will not be made explicit in this note. For example the methods used in the 
percolation version of these arguments [4], see also [5] (with suitable modifications for boundary conditions if 
necessary) can be applied to every (uniform) system known to the authors. 




(1.1) 



th 2 (3J > th0U. 



(1.2) 
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The fact that systems described by Equations (1.1) and (1.2) have a continuous transition is, 
perhaps, not too surprising. Nevertheless, there is only a small handful of models for which a 
continuous transition can be rigorously established. Now it is slightly larger handful. 

2. The Random Cluster Representation. To implement the random cluster representation, 
we rewrite the Hamiltonian in Potts form: 

-m = J2i D (<W; + <W " V6 0iai 6 rtTi ] (2.1) 

with D = 2/3 (J + U) and V = A(3U. For each bond, we may expand as in the usual random 
cluster model, e.g. exp [p8 ai(7j ) = (l + -Rd^^) , etc. with Rp = e D — 1 but when we get to 
the four-spin term, we are forced to acknowledge that exp {—VS ai<Tj S TiTj ^ = (l — ryS aia .S TiT .) 
with — ry = (e~ v — l) a negative number. This could ostensibly lead to negative numbers in 
the graphical representation - and hence no representation with a probabilistic interpretation. 
However the situation is not nearly as bad as it looks since, on each bond, all we are interested in 
is the product of all three terms. Indeed, it is seen that there are there are several mechanisms 
for generating a S aiaj 5 TiTj term and we get 

exp {D (S aiaj + 8 TiTj ) - VS ai(7 .S TiT . } = 1 + Rd^o.oj + RdSt^ + ^VrxW.,^^; (2.2) 

with V\ VD = [R 2 D - r v (l + R D f] ■ Thus, to get started on a decent graphical representation all 
we need is that 9Vr> > 0; i.e. (l — e~ D ) 2 > 1 — e~ v . After some small amount of work, this is 
seen to be precisely the condition (1.2) mentioned in the introduction. 

There is an obvious way to develop a graphical representation. Following closely the derivation 
in [9] for the Potts model, we could define the i?i)5 (Ji(Jj and Rd^t^ terms as single bonds in 
the cr-layer and r-layer respectively while the 9V_D-terms represent double bonds. Then all the 
different types of bonds are treated as separate entities. This was done in [7] and is well suited 
for certain purposes but none that are related to the present work. An alternative approach, 
used in both [8] and [7], is to define bond configurations Q a and Q T for the a and r layers that 
do not care where the bonds came from. I.e. CT (6) = 1 means that b got a single cr-bond or a 
maybe a double bond. In this approach, the resulting random cluster measure is defined by the 
weights 

W(Q a , n T ) = A^ n ' vn ')' B [N(^An T )] 2 c(^) 2 c(n T ) (2 _ 3) 

with A = Rp and B = ^yo/Ro- As usual: N(-) is the number of bonds of the specified type 
while C(Q a ) and C(Q T ) individually count the number of connected components in the cr-layer 
and r-layers respectively. (Of course the latter objects - as well as possible additional constraints 
- are not completely defined until one pays heed to the boundary conditions.) 

It is also convenient to write the weights in the form where each occupied bond gets its own 
factor of A and then an additional factor of B/A for the number of overlapping bonds. Then, 
as is is clear, if B > A, then the bonds in the separate layers attract while if B < A, they repel. 
And notice that the dividing line is determined by the sign of U. What has impeded progress 
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so far is the (straightforward) result concerning the FKG-properties of these measures that was 
proved in [7,8] (see also [10]). For future purposes, we state the result in general terms (c.f. the 
remark following Proposition 2.1) 

Proposition 2.1. Let Q denote a finite graph. For r, s > 1 let / Up C (— ) be the measure on pairs 
of bond configurations (fl a , Cl T ) with weights given by 

b 

where 0, a (b) = or 1 corresponds to the a -bond being vacant or occupied and similarly for 
Q T (b) and P is notation for the various parameters in the weights. Then the sufficient - and the 
necessary - condition for the measure /Up C (— ) to be strong FKG is that > 1 for all b. 

Proof. This is a reasonably straightforward FKG lattice condition to check; see [7,8] (and also 
[10]). □ 

Remark. The definition of the measures - for the benefit of Propositions 2.1 and 2.2 - is a little 
unwieldy. The case in hand is simply A a (b) = A T (b) = A and = B/A. (Thus the theorem 
fails for models of interest in this paper; we will return to this point after we finish apologizing 
for the statement of the proposition.) The primary reason we use "the most general graphs" 
is as follows: to complete the argument that gets us from a Lebowitz inequality to a divergent 
susceptibility we need to consider these measures in the context of sensible graphs - like finite 
pieces of Z 2 - in the presence of boundary conditions e.g. "wired" . On a given finite graph most 
boundary conditions that lead to FKG measures for random cluster problems can be realized as 
free boundary conditions on a related graph. 

As just mentioned: with the usual ordering, the FKG property for the systems described by 
the Hamiltonian in Eq.(2.1) will only hold in case — V > 0. However, there is an alternative 
ordering which is described in [8] that is well suited for the situation at hand. This alternative 
exploits the competitive nature of the interaction between the layers. In particular, for any bond 
6, we can define the local order Q a (b) = 1 > 0, a (b) = and Cl T (b) = > Cl T (b) = 1. Notice that 
with this ordering, there is no relation between the configuration with (Q a (b) = l,O r (6) = 1) 

and (n a (b) = o,n T (b) = 0). 

As obvious as all this may seem in hindsight, it is nevertheless surprising: If, for example, 
U < J, the low temperature states are ordered in both layers: Despite the competition, there is 
ultimately some cooperation. In any case, it is equally not hard to show that for B < A, the 
measures have the FKG property with respect to this alternative partial ordering. 

Proposition 2.2. Let /Up C (— ) be a random cluster measure as described in the statement of 
Proposition 2.1 with < 1 on each bond. Then fip (— ) is FKG with respect to the partial order 
where Q a (b) and 1 — O r (6) are (separate) increasing coordinates. 

Proof. See [8] Lemma 4.1. □ 

The significance of this representation for the present work is exemplified by the following 
(which already appears in [8]): 



4 



L. CHAYES AND K. SHTENGEL 



Corollary. For the spin systems described in Proposition 2.2 let i, j , k and i denote any sites 
of the lattice. Then, for r = s = 2, 

For the other spin systems of this type, the same holds with GiCj replaced by (S ai (Tj — etc. 

Proof. According to the standard translations between spin systems and their random cluster 
representations ([11] or, better yet, [12]) the object (S aiaj — ^)/3n is proportional to the probability 
that the sites i and j belong to the same Q a - cluster. Let us define the various connectivity events: 
{i < — ► j} is the event that i and j are connected in the cr-layer etc. Then we have, for r = s = 2 

{a l a 3 T k T i )f 3 n = ^ j} H {k <—+ (2.5) 

(and similarly for the other values of r and s). However 



ij$fc({i<-^j}n{k<r 




by the FKG property and the desired result is established. □ 

We will use the above inequality - as well as some related inequalities - when the time comes 
to derive our principal results. 

3. The Loop Representation. We now start down an entirely different track namely that 
of loop expansions. The primitive versions of these expansions are well known (e.g. they are 
described in [8]) and, in the context of the present model, actually date back to [6]. However, the 
usual versions will only take us so far and some small modifications will eventually have to be 
inserted. The unfortunate general limitation is that these expansions are heavily tied with the 
Ising nature of the spin variables (or, more precisely the ability to represent the spin variables in 
terms of Ising variables). Thus, for the remainder of this paper, we will confine attention to the 
standard (r = s — 2) Ashkin-Teller models. 

We start with the famous Ising-spin identity: e KO " iCr ? oc l + (th K)OiOj. In any Ising-type system, 
this can be applied to any set that has an Ising-type interaction. But we remark that schemes 
based on this identity usually require the coefficients of the Ising spin terms to be positive. For 
the case at hand we get 

^[Jia^+T^-Ua^nr,] ^ ^ + J^.)^ + J TiTj ){\ - Va^^Tj) = $ itj (3.1) 

With J = th f3 J and U = th (3U. If we multiply these out, the coefficient of the four-spin term 
is seen to be proportional to J 2 — U which we want to be positive. This is, again, the condition 
(1.2). Multiplication of all terms leads to a couple of two-spin terms, a four-spin term and a 
constant (that gets rescaled to one). The result is 



OC (1 + L(7i Oj + hTiTj + KaiO-jTiTj) 



(3.2) 



Lebowitz Inequalities for AT Systems 



5 



with L = [J - U]/[l - J 2 U] and K = [J 2 - U]/[l - J 2 U]. The further development of the 
loop expansion is somewhat obscured by the presence of the four-spin term so to simplify our 
arguments, we will start with the treatment of the case K = - saturation of the condition (1.2). 
(As a matter of fact, this very model was discussed in [13] where, for example, the 2<i-version was 
shown to have a random surface interpretation.) We assume that we are on some finite graph 
that has free boundary conditions. Let us start with the partition function: We get to pick one 
term from for each bond of the lattice. One such choice for each bond leads to a graphical 
configuration (or diagram) which we will call rj: a red bond at if we choose the Lr^Tj-term, 
a blue bond if we pick h<ji<jj and vacant for the "1". When we perform the trace (normalized 
for convenience) the only diagrams that survive are those where an even number - possibly zero 
- of bonds of both types are incident at each vertex. When these constraints are satisfied, the 
weight of the configuration is just HJ 77 '. It is of some use to consider S, a colorless configuration, 
that represents all the configurations rj that have the same bonds regardless of color. Although 
it is tempting proclaim that the weight of H is given 

L |H| x 2 #of loops in the number of loops 

in a configuration is not very well defined - certainly not in any way that makes this formula 
true. (Notwithstanding, we will still call this the loop representation.) In the end, the coefficient 
of lJ~l is just the number of ways that the loop clusters can be consistently colored red and 
blue in such a way that the constraints are satisfied at each vertex. The trivial, but important 
observation is that this counting procedure factors over the individual clusters of H. 

The spin-spin correlation function is not particularly pretty in this representation - it is not 
the probability of anything in the loop measure (i.e. the measure defined by the weights of the 
loop configurations). 3 When we calculate the correlation between o~i and aj the best description 
is to say that there is a numerator and a denominator. The denominator is just the partition 
function - the sum of all the weights described above. The numerator contains similar sorts of 
configurations except that now the two sites i and j have to have an odd number of incident 
(blue) bonds. Since every other site has an even number of bonds, this implies that there is a 
walk from i to j. Following the language of [14] we say there are sources at the sites i and j. 

Let us now consider the four-point function (aiajakCg). Here we have sources at i, j, k and t. 
When we look at the contributing terms, we see that there are two distinct possibilities: All four 
sites are part of the same connected cluster (regardless of color) or the sites pair up forming two 
distinct clusters - and there are three ways of doing this. We denote the sum of the corresponding 
terms by Tij^,e for the case where the four points are tied up in the same cluster, Q[ij\k,e] for % 
and j in a separate cluster from k and t etc. Thus we may write 

Z-r. x (o~iO~jO- k a£)Ti = [T i;jt k,£ + Q[i,j\k,£] + 9[i,k\j,e] + 8[i,e\j,k]\ (3-3) 

where 1Z is the interaction that has J 2 = U. 

Now let us consider a mixed four-point function, say (o~iajTkTe)Ti. Several things are immedi- 
ately clear: First, there are no terms like the last two for this correlation function - we need the 

3 In general the correlation functions can be expressed as the expected value of some function of H — but this 
does not seem to serve any tangible purpose. For certain correlation functions we can do better. For example we 
will show that (<Ji<JjTiTj) is just the probability that i and j belong to the same loop (c.f. Corollary I to Theorem 
3.4). Notice that if U = the former is equal to (aiaj) 2 . This is a disguised version of the famous "switching 
identities" derived in [14]. 
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proper sources to be paired. Second, there will be terms like the first two - a term with all four 
points connected and one with (i & j) and (k & £) separately paired. Thus we may tentatively 
write Z n {ai(jjT k Ti)K = Tij, k ,e + &[i,j\k,£\- Third, it is clear that Q[ij\k,£] = 9[i,j\k,t}- Indeed in 
any colored graph that contributes to Q[ij\ k ,t], one can simply reverse the colors of the cluster of 
(k & £). This shows that 8[i,j\k,e] > Q[i,j\k,e\ and the reverse argument shows Q[ij\k,e\ > Q[i,j\k,e]- 

What is not immediately obvious is that the four-point cluster terms are equal. This is the 
subject of the argument below. 

Lemma 3.1. Let T denote a finite connected graph that contains among others four points k 
and i. All sites in T are evenly coordinated except the above mentioned which are odd- coordinated. 
Let ■ k ^ denote the number of ways that the bonds of this graph can be colored red and blue 
such that an even number of red and blue bonds are incident at each vertex save for i,j,k and 
t which have an odd number of blues. Let jy r ^ k ^ denote the number of similar colorings 
except this time the sites k and t have an odd number of reds. Then 

^b(i,j,k,£) ~ ly b(i,j);r(k,£)- 



Proof. Let p : k — > i denote a self-avoiding path in T with terminal points k and £. Let .A/X ■ k ^ 
be the set of successful colorings of T according to the first set of rules (which is observed to 
be non-empty) and similarly for Af^ jy r ( k iy We will use p to define a map from Af^ j k ^ to 
■^b{% j)-r(k t) which is one-to-one; the map is simply to reverse the color of each bond along the 
path. 

It is clear that the application of this map to an element of A/^ ■ k ^ produces a coloring of 
the desired type: Indeed, at internal sites of the path, the change in the number of blues is ±2 or 
zero with a corresponding change in reds of =f2 or zero. Similarly at the endpoints the number 
of blues changes by ±1 accompanied by a change in the number of reds by =Fl. 

Next, consider two distinct colorings of T in j k iy Either the two colorings have some 
difference on the compliment of p - in which case the target configurations are surely different, 
or they must be different along the path itself - in which case they are still different when their 
paths get reversed. 

Hence this map is indeed one-to-one which establishes NT r /, »>. > NT r ■ u e\- But we can 

1 °(i,3);r(k,t) — b(i,j,k,e) 

construct a similar sort of map from Af^ jy r ( k ^ to Af^ . k ^ which gives us the reverse inequal- 
ity. □ 

As an immediate consequence we get that = Tij jk ,£- This leads to the following 

preliminary result: 

Theorem 3.2. Consider any graph Q with free boundary conditions and interaction 

-1Z = (3 y^(UtCrj + TiTj) - a ViO-jTiTj 

i,3 i,j 

with th a = th j3 . Then 
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Proof. On the basis of what has transpired we have 

r 

(a l a J a k a e )n = (a l a J T k T i ) n + {a^kTjT^n + (a l o-eT j T k )n - 2 l £ k ' . (3.4) 

So we neglect the final term and use the FKG property from Section 2. □ 

Of course the above is of no immediate benefit for the analysis of the susceptibility; in this 
model, there are additional terms when we differentiate the two-point correlation function with 
respect to the temperature. However, the extra ingredients required are not particularly difficult: 

Theorem 3.3. Consider the system described in the statement of Theorem 3.2. Then 

{o-iCFj; [cFkcri + rtrg - a'akO-fTkTiDn < (a i ak)iz(o'jO-e)n + {o-icri)n{o'jO-k)n- 
where, in general, (A; B) = (AB) — (A)(B) and a' = da/d/3. 

Proof. We first remark that all we really require is that the coefficient of the four-spin term is 
not less than —1. We observe that in the present case, a' = 2th/3/[l + th 2 /3] < 1. In light of 
Theorem 3.2, all that we need to show is that (o~iO-j; \jkJt — ct ' OkO(TkT(^Ti is not positive. Using 
a' < 1, we start with (1 — a')(aiO~j;TkTg) < by the FKG inequality. Now what is left is to 
show that (o-iO-j] [rkT£ — OkOUkT^n < 0. We go back to the FK representation and consider the 
untruncated quantity (<Tja jTkTi(l — o~kcrt))fc. It is clear that this vanishes unless k is connected 
to £ by r-bonds. But this also vanished unless k is not connected to £ by cr-bonds. If both of 
these are satisfied, and the site i is a-connected to j, we get one. There are other possibilities 
for getting something non-zero but these all come out negative. We thus have 

(cno-jinn - a k a e T k T e ]) n </u| C ({«^j'}n [{k <— > £} n {k (3.5) 

where {k <—x^> £} is notation for the the compliment of the event {k < — ► £}. The second event 

a cr 

on the right hand side (the one in the square brackets) is clearly decreasing while the first one is 
increasing. Thus 

{<Ti<Tj[TkTt - o- k aiTkT £ ])n < Vn C ({* <~^> j}) /4?({ k ^ £} n {k £]). (3.6) 

The second factor on the right hand side of Eq.(3.6) is seen to equal (tv-t^I — o-kO~i\)n and the 
desired version of the Lebowitz inequality is established. □ 

The corollary to the above should be that the susceptibility diverges. However before such a 
claim can be made, it has to be verified that the susceptibility is actually infinite somewhere. 
Unfortunately with the above interaction, this can only be established for the model on the 
square lattice; here it seems likely that here the critical value is in fact (3 = oo (c.f. the discussion 
in [13]). 
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Corollary. For the model on Z 2 with the interaction as described in the statements of Theorem 
3.2, the susceptibility diverges continuously at some value - possibly infinite - of the parameter 
P- 

Proof. Infinite susceptibility, at (3 = oo, was established for this model in [13] using duality, 
the rest follows from the bound on dX/d/3 that results from the inequality in the statement of 
Theorem 3.3 □ 

As emphasized by the above limitations, it is clearly desirable to extend this set of results past 
the K = restrictions. Going back to Eq.(3.2) - with K > - the straightforward procedure 
would be to include one new type of bond, namely a red-blue (double) bond that has weight IK. 
As indicated, these bonds have double occupancy but all the other constraints are pretty much 
the same as before. Although this is the standard way of doing these loop expansions, it is clear 
that our analysis will run into trouble. 

First, on the aesthetic level, there is now no worthwhile expansion for the model in terms of 
colorless graphs. Every double bond forces two colors through that bond. This means that under 
certain circumstances, the color of one bond in a cluster could determine the coloring scheme of 
the entire cluster - even though the cluster itself could have many "loops" . On a more disturbing 
note: Although an expression like Eq.(3.3) is still valid for the usual four-point function, and 
similarly for the mixed four-point functions, there is no obvious correspondence between the 
various T-terms. (We remind the reader that these are the terms corresponding to the diagrams 
where all four points are in the same cluster.) Indeed it is clear, after a moments reflection, that 
each of these T's contain diagrams that are conspicuously absent in the other. (The relevant 
0-terms do coincide but this is small consolation.) The resolution of all of these difficulties is, 
perhaps, obvious in hindsight: We will borrow from the one. 

What we actually do is as follows: We write 1 = (1 — K) + K (noting that IK < 1) so that 

3>ij OC ([1 — IK] + JOiOj + §TiTj + KaiO-jTiTj + IK) 

OC (1 + J(TiO-j + JTiTj + KaiO-jTiTj + K) (3.7) 

with K = K/[l — K] etc. The "new" terms corresponding to the selection of the K may, perhaps, 
be called dud-bonds. In the colored graphs these have no special significance - they do not 
change the color parity of any vertex - but in any colorless representation, they have the same 
appearance, not to mention the same weight, as the red-blue type of double bond. 

This is the bottom line, anything else may be regarded as window dressing. Notwithstanding, 
an interpretation that is pleasing (at least to the authors) is as follows: Each edge of the lattice 
can potentially be split into two edges, an "upper" edge and a "lower" edge. Let us write the K 
terms in Eq.(3.7) as K = ^K rr + ^K bb and KaiajTiTj = ^K. rh aiajTiTj + ^K. hr aiajTiTj with (of 
course) K = K bb = K rr = K rb = K br . Now, in a double bond situation (which automatically 
implements the "bond splitting" option) one of these four terms is selected. The four choices 
are identified as double blue, double red, upper-red/lower-blue and upper-blue/lower-red i.e. 
the two edges of the split actually colored. The constraints are exactly as they were before: an 
even number of colors of each type at each vertex. With this interpretation in mind, it is clear 
that we once again have a bonified - and useful - colorless representation. Explicitly, for a given 
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configuration, the graph becomes a graph with single edges and split edges. The weight is [a 
factor of J for each (occupied) single bond and a factor of for each (occupied) double] x [the 
number of ways that the configuration can be consistently colored]. Notice that in any coloring 
in which a particular double ends up red/blue, there is another coloring - regarded as distinct 
- with all other bonds the same but the particular bond a blue/red. Similarly for exchanging 
a red/red with a blue/blue (but notice that the mono-colored and di-colored doubles cannot be 
exchanged without forcing other changes in the configuration.) Thus, when all is said and done, 
we are back to the situation we had at K = - albeit on the peculiar graph defined by the splits. 
As a consequence obtain, rather inexpensively, 

Theorem 3.4. Consider a finite graph with the Hamiltonian 

with th 2 Ji j > th(3Uij > on every edge Then 

(o io jo ko -e) pn < {o~iO j) pn{o 'ko e) pn + (ciCfc) pn(o 'j a e) pn + (o~iO-i)pn(o~jO-k}pn- 

and 

(o-iO-j-, \jkJt - o ko eTkTe]) pn < 

Proof. The second inequality is a direct consequence of the FKG property (that was established 
in Proposition 2.2) by following the argument in Theorem 3.3. Thus we are left with the extension 
of Theorem 3.2 to non-zero K. Let us start by repeating Eq.(3.3) and the similar equation for 
the mixed four-point function without notational changes. Obviously Q[ij\k,i] still equals Q[ij\k,t] 
by blatant color-reversal symmetry. Of more significance, we claim that with the dud-bond 
representation, (and the ensuing interpretation) the extension of Lemma 3.1 is immediate. In 
particular the argument that was used to show Tij^i = Tij } k,i (whereby some path connecting 
k to £ gets color reversed) still goes through - the formulation was for a general underlying graph 
and thus works for the graph defined with the splits. Given the identities V^j^i = ^i,j,k,e and 
8[i,j\k,e] = 9[i,j\k,e] ( as wen as the FKG property) the desired result is manifest. □ 

As a bonus, we obtain a direct relationship between a correlation function and a loop proba- 
bility: 

Corollary I. Consider the system as described in the statement of Theorem 3.4- Let Cij = 
Cij(PTi) denote the probability, in the loop representation, that the sites i and j belong to the 
same loop cluster. Then 

Proof. The numerator of the right-hand side requires a blue walk and a red walk between the 
sites % and j. Applying the color switching techniques described in Lemma 3.1 in the context of 
the dud-bond representation - as was done in Theorem 3.4 it is seen that each such term is in 
one-to-one correspondence with a term contributing to Cij and vice versa. □ 
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Corollary II. Consider a periodic graph on which there is Hamiltonian TL (also periodic) of the 
type described above with Uij = (3~ 1 Aij((3) where the Aij((3) are smooth functions satisfying 
th Aij(f3) < th 2 (3Jij as well as ■(/?) < Jij and lim^oo (3~ 1 A i j(f3) < J itj . Here (3 < (3 < oo 
with Pq small enough to ensure that at f3 = (3q, the susceptibility is finite. Then there is a f3 c 
along this trajectory where the susceptibility diverges continuously; in particular at least as fast 
as a constant times {(3 C — 

Remark. One can presume that under the stated conditions every ferromagnetic transition in the 
region U > is covered. Although this may not seem evident in the physical variables used above, 
it is quite transparent in the A—B variables that were featured in Section 2. For example, in the 
homogeneous case with J = 1 and j3U = a, let us set B = XA with A < 1 and allow A to range in 
[0, oo). Then we find that r v = (1-X)[A/(1 + A)] 2 and, in general, a' = r v /[l-r v (l+A)}. Thus 
a' < 1 because the worst case is the one with the biggest value of ry (namely A = 0) which is the 
borderline case discussed in Theorem 3.3. Now starting at A = there is exponential clustering 
for all A small enough. But for A large and A > the graphical representation discussed in 
Section 2 has percolation (in both layers) which implies magnetization. (These facts can be 
verified by simple-minded comparisons to independent percolation.) Thus, along these curves, 
there is always a point where the susceptibility is infinite. So, for all intents and purposes, this 
region has a phase boundary separating (3 = from (3 = oo. For negative values of B (i.e. A < 0) 
the ground states are not ferromagnetic and it seems highly unlikely that ferromagnetism could 
somehow be generated by entropic effects. 

Proof. The above hypotheses are sufficient to derive bounds of the form dX/d/3 < [const. \X 2 
(suitably regularized in finite volume if necessary). Furthermore, they are sufficient to ensure 
a ferromagnetic transition at finite (3 (c.f. the above remark) and hence X is actually infinite 
somewhere. The combination of these ingredients implies the divergence of the susceptibility 
with the stated bound. □ 
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